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In this work, we present a study of a purely kinetic k-essence model, characterized basically by 
a parameter a in presence of a bulk dissipative term, whose relationship between viscous pressure 
n and energy density p of the background follows a polytropic type law 11 oc where A, in 

principle, is a parameter without restrictions. Analytical solutions for the energy density of the 
k-essence held are found in two specihc cases: A = 1/2 and A = (1 — a)/2a, and then we show 
that these solutions posses the same functional form than the non-viscous counterpart. Finally, 
both approach are contrasted with observational data from type la supernova, and the most recent 
Hubble parameter measurements, and therefore, the best values for the parameters of the theory 
are founds. 

PACS numbers: 04.20.Fy, 04.20.Jb, 04.40.Nr, 04.70.Bw 


I. INTRODUCTION 


At present, the scientific community dedicated to the 
study of the universe have deep and intriguing questions 
unanswered. One of the most fascinating corresponds to 
what we know as dark energy (DE) [lH5[, a component 
designed to explain the current acceleration in the ex¬ 
pansion of the universe. In its simplest form, this can be 
described by a perfect fluid with constant energy density, 
which leads to the useful A - cold dark matter (ACDM) 
model, the simplest model that fits a varied set of ob¬ 
servational data. However, this model has a high de¬ 
pendence to initial conditions that makes it unnatural 
in many ways. For example, the current value for Da 
and VloM are of the same order of magnitude, a fact 
highly improbable, because the dark matter (DM) con¬ 
tribution decreases with a~^, with a{t) the scale factor, 
meanwhile the cosmological constant remains constant. 
This problem in particular is known as the cosmic co¬ 
incidence problem. It is for this reason that many of 
the most sophisticated experiments and instruments have 
been put in place; as the Dark Energy Survey (DES) 
a, the Baryon Oscillation Spectroscopic Survey (BOSS) 
7|, and the upcoming Large Synoptic Survey Telescope 
(LSST) 0 to mention some, all of them trying to find 
new insights into the nature of dark energy. In this con¬ 
text, the most natural way to understand the acceleration 
of the universe, is to assume the existence of a dynamical 
cosmological constant, or a theoretical model with a dy¬ 
namical equation of state parameter {p/p = w{z)). The 
source of this dynamical dark energy could be both, a new 
field component filling the universe, as a quintessence 
scalar field [m, or it can be produced by modifying 
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gravity Ilw22ll . In this work, the so-called k-essence 
model 0] |24| is used, which is a type of dynamical cos¬ 
mological constant model, but where the source of its 
dynamics comes from a non trivial kinetic term, as op¬ 
posite to the case of a typical quintessence model where 
the source is a different scalar field potential, and then 
put it into the test with current observational data from 
both, type la supernovae [1^ , and the most update Hub¬ 
ble parameter measurements [2^. 

Besides, if we focus, for example, in the dark sector as 
a whole, it has been proved that the division of this sec¬ 
tor into DM and DE is merely conventional since exist a 
degeneracy between both components, resulting from the 
fact that gravity only measures the total energy tensor 
M (see also (281 - 0 ). So, in the lack of a well confirmed 
detection (nongravitational) of the DM only the overall 
properties of the dark sector can be inferred from cosmo¬ 
logical data, at the background and perturbative level. 
This results has driven the research to explore alternative 
models which consider a single fluid that behaves both 
as DE and DM, the called unified DM models (UDM). 
So this fluid must drive both the accelerated expansion 
of the Universe at late times and the formation of struc¬ 
tures (see 0 for a review of these models). Of course, 
a small speed of sound should be an essential character¬ 
istic of a viable unified model in order to do not impede 
the structure formation and to have a ISW effect signal 
compatible with CMB observation (36l - l4ll| . 

In this present work we will consider UDM models de¬ 
rived in the framework of k-essence fields, common in 
effective field theories arising from string theory and in 
particular in D-branes models [421 - 10 . This generaliza¬ 
tion of the canonical scalar fields models can give rise to 
new dynamics not possible in quintessence. In the con¬ 
text of cosmology, k-essence was first studied as a model 
for inflation (k-inflation) [13] . K-essence models has also 
addressed the problems of a dynamical DE [ 13 , | 13 | and 
the coincidence problem [13, [5l|. For example, a partic¬ 
ular case is the Generalized Chaplygin gas (GCG) which 
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ap pea rs as the simplest tachyon field model, introducing 
in [^, with a constant potential. Moreover, k fields leads 
to a new Chaplygin gases. Within the models investi¬ 
gated in order to unify DE and DM are the 
and those known as purely kinetic models l62l|. The 
unification of DE, DM and inflation has been addressed 


Another issue that emerges from the cosmological data 
is that the exotic behavior of the universal fluid can be 
characterized by a negative pressure and usually repre¬ 
sented by the equation of state w = pjp, where w lies 
very close to —1, most probably being below —1. For 
example, the last Planck results give w = — l.lSto ig 
and w — —1.090.17 (95%CL) by using CMB combined 
with BAO and Union2.1 data [ 63 , respectively, for a con¬ 
stant w model. In combination with SNLS3 data and 
Hq measurement, the EoS for this dark component are 
w = — 1 . 13 f'^Q;i 4 and w = — 1 . 24 ^Q; 4 g {2aCL), respec¬ 
tively. The possibility of ic < — 1 is favored at the 2cr 
level. These results are indicating that a phantom be¬ 
havior of the dark energy component can not ruled out 
from current cosmological data. 

As it was pointed out in [o^ dark energy with a con¬ 
stant EoS ra < — 1 leads to uncommon cosmological sce¬ 
narios. First of all, there is a violation of the dominant 
energy condition (DEC), since p + p < 0. The energy 
density grows up to infinity in a finite time, which leads 
to a big rip, characterized by a scale factor blowing up in 
this finite time. Nevertheless, sudden future singularities 
are not necessarily produced by a fluids violating DEC. 
Solutions which develop a big rip singularity at a finite 
time without violate the strong-energy conditions p > 0 
and p -f 3p > 0 were found in [^ . 1^ . Studies of uni¬ 
fied dark matter models, which are generalizations of the 
Chaplygin gas, presents EoS w < — 1 but without a big 
rip type solution in 


Another mechanism that allows a violation of DEC 
is the existence of dissipation within the cosmic flu¬ 
ids [zainl. In the case of isotropic and homogeneous cos¬ 
mologies, any dissipation process in a FRW cosmology is 
scalar, and therefore may be modeled as a bulk viscosity 
within a thermodynamical approach. The bulk viscos¬ 
ity introduces dissipation by only redefining the effective 
pressure, p^^^ , according to p^j^ = p -I- B = p — 3((i7 , 
where 11 is the bulk viscous pressure, C is the bulk vis¬ 
cosity coefficient and H is the Hubble parameter, and 
c = 8 ttG = 1 (as in all the work). Since the equation of 
energy balance is p -I- 3H{p -|- p -|- H) = 0, the violation 
of DEC, i.e., p-|-p-|-n < 0 implies an increasing energy 
density of the fluid that fills the universe, for a positive 
bulk viscosity coefficient. The condition C > 0 guaranties 
a positive entropy production and, in consequence, no vi¬ 
olation of the second law of the thermodynamics . 

Some investigations have considered that the viscous 
pressure can drives the present acceleration of the Uni¬ 
verse, so it can be used to eliminate the dark energy 
component and to formulate unified dark matter model 
with viscous pressure. In [zl,[z3, for example, cosmolog¬ 


ical models where the only component is a pressureless 
fluid with a variable and constant bulk viscosity was con¬ 
fronted with the observational data. Nevertheless, the 
bulk viscosity induces a large time variation of the grav¬ 
itational potential at late times which leads to inconsis¬ 
tencies with the integrated Sachs-Wolfe (ISW) effect in 
such model In order to overcome this problem, 

Velten & Schwarz proposed a model with a viscous 
cold dark matter and a cosmological constant, which acts 
driving the accelerated expansion of the Universe. Our 
aim in this work is to investigate UDM models derived in 
the framework of k-essence fields which can also present 
dissipative effects. 

Usually k-essence is defined as a quintessence, scalar 
field (j) with a non-canonical kinetic energy associated 
with a Lagrangian £ = —V{(j))F{X). In the subse¬ 
quent calculations, we shall restrict ourselves to the sim¬ 
ple k-essence models for which the potential V = Vq = 
constant. We also assume that Vb = 1 without any 
loss of generality. One reason for studying k-essence 
is that it is possible to construct a particularly inter¬ 
esting class of such models in which the k-essence en¬ 
ergy density tracks the radiation energy density during 
the radiation-dominated era, but then evolves toward 
a constant-density dark energy component during the 
matter-dominated era. Such a behaviour can to a cer¬ 
tain degree solve the coincidence problem. 

We investigate a dark energy model described by 
an effective minimally coupled scalar field with a non- 
canonical kinetic term. If for a moment we neglect the 
part of the Lagrangian containing ordinary matter, the 
general action for a k-essence field (j) minimally coupled 
to gravity is 


S = Sg + = - J , (1) 

where F((j)^X) is an arbitrary function of (j) that rep¬ 
resents the k-essence action and X = (j) is the 

kinetic term. We now restrict ourselves to the subclass 
of kinetic k-essence, with an action independent of (j) 

S^ = -j d^x^F{X). ( 2 ) 

Unless otherwise stated, we consider (j) to be smooth on 
scales of interest so that X = > 0. The energy- 

momentum tensor of the k-essence is obtained by varying 
the action ([2]) with respect to the metric, yielding 

= Fxd^fjiducj) - gfj.uF, (3) 

where the subscript X denotes differentiation with re¬ 
spect to X. Identifying as the energy-momentum 
tensor of a perfect fluid, we have the k-essence energy 
density p and pressure p 


p = F-2XFx, (4) 

p = -F (5) 
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Throughout this paper, we will assume that the energy 
density is positive so that F — 2XFx > 0. The equation 
of state for the k-essence fluid can be written as p = 
w^p = (70 - l)p with F > 0, 

= = l = 2XF„-F - ® 


de Sitter stage, inducing a modified Chaplygin gas. The 
explicit equation of state can be obtained from Eqs. (HD 
and ([5]) 

p= (n-l)p-naop^, (15) 

where the parameter n is a function of the constant a, 
given by 


II. THE K-ESSENCE MODEL WITH 
DISSIPATION 


2a 

2a-I 


(16) 


The Friedman-Lemaitre-Robertson-Walker (FLRW) 
metric for an homogeneous and isotropic flat universe 
is given by 

ds^ = -df + a{tf [dr^ + {d9^ + sin^9 dif^)] , (7) 


Obviously, the range of this parameter is 1 > n > 0, 
if —00 <a<0;0>n> —oo, if 0 < a < 1/2; and 
oo>n>l, ifl/ 2 <a<oo. 

Of course, the speed of sound is affected by the viscous 
pressure, which becomes 


where a{t) is the scale factor and t represents the cosmic 
time. In the framework of the first order thermodynamic 
theory of Eckart the field equations in the presence 
of bulk viscous stresses yield 





P 
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( 8 ) 
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where is the speed of sound in the purely k-essence 
background [l^, given by 


V 


2 

0 


(n 




(18) 


- — H + — —- (p -I- 3p , (9) 

where the effective pressure is given by 

Pe//=P + n, (10) 

and 

n = -3i7C, ( 11 ) 

is the bulk dissipative pressure and the viscosity. In 
what follows we will assume a power law dependence for 
the viscosity in terms of the the density 

C = CoP^ (12) 

where Co is a positive semi-definite constant with dimen¬ 
sion T~^, and A may take any value. For 

example, the most common values are A = 1/2, i.e., 
C oc [soHs^l and A = 1, i.e., C oc p [11,11^. These 
values were chosen because leads to well known analytic 
solutions. Therefore, the conservation equation for the 
fluid can be written as 


From Eqs. (El-ITll), together with the EoS we ob¬ 
tain the evolution equation for id in terms of the redshift, 

-ao^+aiii + a2ii''-^+a3id^-^ = 0, (19) 

ax 

where x = ln(l + z)., and the coefficients are given by 

Oo = 2, oi = 3n, 02 = -3“as = -3'^’'‘^Co, (20) 
whereas the exponents reads 

P = 2(^^^^=b /3 = 2A-fl. (21) 

As a first observation, we note that there are two spe¬ 
cial values that yields to well know equation without vis¬ 
cosity [i^: /3 = 2 (A = 1/2) and /3 = p (A = ^ 1/2). 

These values leads to a single equation which posses a 
generic structure for its quadrature given by 

(ITT 

— = Ai ii + A2 id"-^ = Ai (id + y id''-!) , (22) 


p + 3H{p + p + n) =0. (13) 

In this work we consider the following function F for the 
k-essence field 

F{X) = - 2aaQ^/X], (14) 

2a — i 

where a and ag are two real constants. This generating 
function exhibits a transition from a power law phase to a 


where y = A^jA-y^ and the new coefficients are given in 
terms of the above by the following expressions 


Ai 


ai + 03 
Oo 


^(n- TsCo), 


Ai 

y = 


flo 

ai -I- 03 


Q 

3 It ofQ n 

2 ’ 

_ gp n 

■ 3i(V3Co-n)’ 


(23) 

(24) 

(25) 
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for A = 1/2 (the model A), whereas 
3 


Ar = ^ = 


Ao = 


gq 1 — 2A 
02 + 03 3 ^+^ / 2ao 


y = 


oq 2 

02 + 03 _ 3^“^ 

Gl 2 


(26) 

(y3Co(2A-l)-2ao) ,(28) 


for A 7 ^ 1/2 (the model B). So, a direct integration of Eq. 
(HH) leads to 


Hiz) = Ho 


(1 + z)^-A3 
1 — A3 


(29) 


and therefore, the energy density is given by 


p{z) = iHl 


(1 + z) - A 3 
1 — A3 


where we have defined 


A 3 


7^ 

i + n' 



(30) 


(31) 


We note that the generic expressi on (1501) (or Eq. (HID) 
has the form found by Chimento [52l |. and obviously, 
these case is entirely recuperated by making ^0 0 and 

A —>■ 0. A second observation is that, in the case A 7 ^: 1/2 
and by using Eqs. (EU and ( 1 ^ . the expression (Hi) 
takes the form 


H{z) = Ho 


(1 + z )3 - A3 

I-A 3 


1-2A 


(32) 


Finally, there is a future singular value of the redshift, say 
Zs^ for which the Hubble function takes its zero value: 

2 , = Af^ - 1. (33) 


We are restricted to the realistic values for the future 
singularity, so we expect that —1 < < 0. Thus, this 

condition impose that y > 0, which implies that n < 
V^Co if A = 1 / 2 , and A > 1/2 + ao/('/3Co) if « = (2A + 
1)“^. In this context, notice that A = 1 (i. e. a = 1/3) 
leads to the condition 2 ao < a/SCo- 

This kind of future singularity correspond to a novel 
type, because although both the Hubble parameter (1^ 
and the energy density dsni) vanish at this redshift. 


III. OBSERVATIONAL CONSTRAINTS 

In this section we use observational data to put some 
constraints in the free parameters of the models. We use 
type la supernova data, specifically the Union 2 data set 


[ 2 ^, and the most recent Hubble parameter H{z) mea¬ 
surements compiled in [s^, consisting in 28 data points 
expanding a range in redshift 0.015 < z < 2.3. 

The comoving distance from the observer to redshift 
z, in a flat universe, is given by 


r 




dz' 

Wv 


(34) 


where E{z) = H{z)/Ho. The SNIa data give the lu¬ 
minosity distance dL{z) = (1 -I- z)r{z). Notice that the 
procedure we follow differ from those used by Bandy- 
opadhyay et al. [s^. In this work the authors define an 
intermediate parametrization for the luminosity distance 
as a function of two parameters a, /3 which after the fit¬ 
ting is related to the physical parameters of the model. 
Here we constrain directly the physical parameters of the 
model. 

We fit the SNIa with the cosmological model by mini¬ 
mizing the value defined by 


2 _ [y{zt) - yobs{Zi)Y 

XSNIa — 2 ^ 
i=l 


(35) 


where p{z) = 5 logio[dL(z)/Mpc] -|- 25 is the theoretical 
value of the distance modulus, and pobs is the correspond¬ 
ing observed one. 

From (1301) we can write down explicitly 


E{z) 


- 2Ai 

(1 -I- Z)~i^ — A 3 
1 — A 3 


tx 

2 


(36) 


This form of the solution enable us to test both models 
at the same time by reinterpreting the constants values. 
The best fit values using both SNIa and H{z) data leads 
to a xYed — 0-96, and Ai = I.50±0.15, rj = —0.29±0.19, 
and A 3 = —2.9 ± 0.5. 

For the case A = 1/2, the free parameters are three: 
Al, the parameter that changes with the model, rj, which 
is defined in (EH), and A 3 , defined in m- Straight¬ 
forward calculations lead to n = 0.87 ± 0.07, Co = 
-0.075 ±0.069, and ao = (2.96 ±2.6) x 10“'*. Note that 
for the case A = 1/2 and since we have taken G = I/ 87 r 
and c = I, it is straightforward to see that parameter Co 
is dimensionless. 

Since the exponent in (1551) reduces to = 3, in the 
case A 7 ^: 1/2 the free parameters reduce to A 3 and rj. 
For this reason, is not possible to invert the equations 
completely, because this model is described by three pa¬ 
rameters, Co,Q:o and A. In fact, from the best fit, we 
can write down directly the value for A = —0.65 ± 0.08. 
The other two parameters are tightly related through the 
relation 

«o = -^3^±2V3(2A-1)Co. (37) 

Because the best value for parameter y is large compared 
with the second term in the right hand side (for reason¬ 
able positive values of Co), the value for ao is largely 
better constrained than Co- 
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In order to make manifest the quality of the fit of our 
models, in Figure © we show the theoretical curves of 
best fit for each model together with the observational 
data of H{z). There we show the 28 data points mea¬ 
surements of the Hubble parameter together with the 
best theoretical fit. We have to notice that although the 
lines does not seems to follow the observational points 
very well, this is because the best fit model was computed 
using both SNIa data and H{z) measurements, and the 
first data set statistically weighs more than the second 
one, just because of the number of data in each case. We 
also display in figure ([5]) the confidence level contours for 
the parameters 77 and at one and two a, and in figure 
([3|) the confidence contours for the model B parameters. 



FIG. 1: Using the values of the best fit for each model, here 
we display the theoretical curve of each model along the ob¬ 
servational data for H{z). The continuous line is model A, 
and the dashed line describe model B. It should be noted that 
the values of the best fit was obtained using both measure¬ 
ments of H{z) and supernovas. We have adopted h = 0.673 
from the Planck Collaboration | 88 ||. 

In both cases, because the analysis was performed 
without imposing external priors on the parameters, we 
found a preference for nearly zero to negative values for 
the viscosity constant ^o- 

Despite the strange results - a negative value for the 
viscosity constant - after put in tension our solutions 
with the data, we have confident that such a analysis 
can be done in the first place for any other analytical 
solution that can be obtained in the future. Of course, 
we do not expect to find that just our special (analyt¬ 
ical) solutions be the best fit to the data immediately. 
Cosmology has entered into the era of precision cosmol¬ 
ogy, and with it, the possibility to rule out effectively a 
particular cosmological model. 

IV. CONCLUSIONS 

We have analysed the general relations for a model of 
k-essence generated by the function F{X) = — 


2aaf)\/X] ( proposed by Chimento [l^), when a dissipa¬ 
tive pressure H oc is included. We found a family 

of analytical solutions in two special cases: A = 1/2 and 
A = (I — a)l2a (with a ^ 1/2), which coming from a 
similar differential equations and posses the same struc¬ 
ture that the non-viscous case (compare, for example, 
Eq. (69) in reference with Eq. (pTill ). 

Also, a quick observation of Eq. (HB shows that, de¬ 
pending on the value of A, the speed of sound may be 
greater (A < —1/2), equal (A = —1/2) or less (A > —1/2) 
than the speed of sound without viscosity. Obviously, a 
well behaved fluid requires A > —1/2, which corresponds 
to a consistency relation for A. 

As a light of observational data, we confront both an¬ 
alytical solution with measurements of i?(z) and super¬ 
novas. The best fit yields the following values for the pa¬ 
rameters: xled — 0-96, Ai = 1.50±0.15, p = —0.29±0.19, 
and A 3 = —2.9 ± 0.5. Therefore, we obtain for the 
model A that n = 0.87 ± 0.07, Co = —0.075 ± 0.069, 
and ao = (2.96 ± 2.6) x 10“^, while for the model B we 
obtain that A = —0.65 ± 0.08, and the other two param¬ 
eters are tightly related through Eq. [37l So, both model 
present a controversy with the physical meaning (Co < 0 
in the model A and A < —1/2 in the model B). 

This work can be improved in many ways. On one 
hand, we can attempt an alternative way to obtain an¬ 
alytical solutions, a possibility we are already studying 
using a novel techn ique pr oposed to solve complex dif¬ 
ferential equations [901 - 19^ . In this case, we have the 
possibility to consider A a free parameter, enhancing the 
parameter space to find a best fit with data. 

Certainly a more realistic model would also be in¬ 
teresting to study. In this work we have considered a 
UDM model assuming nothing else but a k-essence field 
is present. We can add explicitly a dark matter term 
and/or a radiation component. We are interested in test¬ 
ing if adding these terms would alleviate our concerns 
about the sign of the viscosity coefficient. 
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FIG. 2: Here we display the 68.27% and 95.45% confidence 
regions for the parameters Ai, As and rj for model A. 
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FIG. 3: Here we display the 68.27% and 95.45% confidence 
regions for the parameters A 3 and 77 for model B. 



